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Abstract 

In this paper we discuss the construction and the dynamics of vortex-like 
topological spin excitations in the Schwinger-boson description of Heisen- 
berg antiferromagnets in two dimensions. The topological spin excitations are 
Dirac fermions (with gap) when spin value S is a half-integer. Experimental 

and theoretical implications of these excitations are being investigated. 
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Recently there are strong interests in the study of Heisenberg antiferromagnets on two 
dimensional square lattice. Early theoretical analysis showed that the low energy, long length 
scale behaviour of the model should be described by the renormalized classical regime of the 
(2 + l)d 0(3) non-linear-sigma-model (NLaM) jl],0. The analysis was believed to be valid 
as long as the length scale under investigation is much larger than lattice spacing. More 
recently, in a series of experiments |§ and Monte Carlo simulations J|-f|, it was discovered 
that the length scale at which the renormalized classical cx-model description becomes valid 
is surprisingly long (L > 200 sites for spin S = 1/2) [|]|§. Quantum Monte Carlo study of 
the low energy spectrum of the model on finite size lattices also suggested that the spectrum 
in the S = 1/2 case disagrees rather strongly with prediction of NLoM even when the size of 
the system N is not too small (N < 32) 0. These anomalous findings suggest that contrary 
to usual believe, there may exist an intermediate energy/length scale where the behaviours 
of the system deviates strongly from NLaM. 

On the other hand, it has been shown by Read and Chakraborty || that fermionic spin 
excitations exist as topological excitations in the short-ranged RVB phase of 2D S = 1/2 
Heisenberg model, and their analysis has been generalized by the author || to the case 
(for arbitrary spin S) when the system is described by Schwinger-boson mean-field the- 
ory (SBMFT) [0. These topological excitations are finite energy excitations of 2D Heisen- 
berg antiferromagnets that do not appear in conventional NLaM description. However, the 
dynamics of the topological spin excitations were not investigated in those works P|| . 

The purpose of the present paper is to investigate, starting from SBMFT, the dynamics of 
the topological spin excitations constructed in Ref. || under some very general assumptions. 
In particular, we shall show that the existence of these topological excitations is consistent 
with the recently observed anomalous behaviours in the Heisenberg model. We shall also 
point out, in the case of S — 1/2 Heisenberg antiferromagnet, an interesting connection 



between the topological spin excitations and the spin excitations in the flux-phase |Tl| of 
fermionic RVB mean-field theory. 

To start with we review briefly the construction of topological spin excitations in SBMFT 



of 2D Heisenberg model. Following Ref. |{J we assume that the low energy physics of general 
spin-5 1 2D Heisenberg model is described by the Schwinger-boson mean-field Hamiltonian 

Hmf = ~JJ2 (&*(Z^Z n - ZnZtf) + B.C. - |A| 2 ) + J2 KZiaZ ia - 25), (1) 

i,j ier 

where A(A*) and A are mean-field parameters determined by the mean-field equations 
A =< (Z^Zji — ZiiZjj) > and < Z^Z^ > + < Z^Z^ >= 2S. Notice that in two di- 
mensions, long-ranged antiferromagnetic order exists at zero temperature corresponding to 
bose-condensation < Z >^ in SBMFT. We shall consider finite temperature T ^ and 
< Z >= in the following. The effect of bose condensation will be addressed at the end of 
the paper. Notice that SBMFT offers a fairly accurate description of the low energy physics 
of Heisenberg model at two dimensions flCfl , thus justifying the starting point of our theory. 



To look for topological excitations in SBMFT, we notice that the structure of the mean- 
field theory resembles very much the BCS theory for superconductivity, except that the spin 
pairs of bosons replace the electron (fermion) Cooper pairs in BCS theory. The resemblance 
of the two theories leads us to study vortex excitations in SBMFT, since vortices are stable 
topological excitations in BCS theory at two dimensions. In BCS theory, a vortex located at 
r = is a solution of the BCS mean-field equation where the order parameter ABcs{r) has 
a form Abcs^) = f(r)e ld in polar coordinate, where f(r) is real and positive. To minimize 
energy, a magnetix flux of 7r-flux quanta is trapped in the vortex core. The vortex solution in 
SBMFT has the same structure, except that the BCS order parameter A BC s is replaced by 
the Schwinger boson order parameter Ay and the vector potential A does not represent the 
physical magnetic field, but is a fictitious gauge field arising from phase fluctuations of order 
parameter Ay's The existence and stability of the vortex solution was demonstrated 



in Ref. || in an effective Ginsburg-Landau theory. These vortex solutions are bosonic, S = 
topological excitations in SBMFT ||. 

To construct topological spin excitations we note that like vortices in superconductors 
where electronic bound states often exist inside the vortex core, boson bound states may 
exist inside vortices in SBMFT. In particular, we argued in Ref. 0] that for vortex centered 



at a lattice site, a bound state of 2S bosons must be formed at the vortex center because of 
the constraint that there is always 25* bosons per site in the Heisenberg model. In particular, 
because of statistics transmutation associated with binding quantum particles to flux-tube 



of 7r-ffux quantum in two dimensions |13| , the resulting excitation is a spin-S 1 fermion when 
2S is odd 0. 

We shall now supply the mathematical details. An important difference between SBMFT 
and BCS theory is that the mean-field Hamiltonian (P breaks translational symmetry of 
the Heisenberg model by one lattice site. Correspondingly there exists two lattice sites per 
unit cell in SBMFT and the fluctuations of the order parameter A^'s are described by two 
amplitude and two phase (uniform and staggered) fields in the continuum limit [0,0], i.e. 
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A,,, l± , = - ± -) + g ± ,(* ± -) J e L J ' 2 'J, (2) 

where q- u = —q u , A s _ v = —A s v are 'staggered' components of the amplitude and phase 
fluctuations of A, respectively. <fi and j x A.dx' are the corresponding 'uniform' components. 
The effective Ginsburg-Landau action for the continuum field variables is derived in Ref. 
H(see also Ref. |12[)- We obtain to order 0(m°) (m is the mass gap for spinwave excitations 



in SBMFT which is very small at low temperature (m ~ Je JS ^ T )), S e ff = S u + S s , where 
S u ~ jdr j d 2 x ((a - 4(S + ±))<f> + 2 + (S + i)0(3 M - , (3) 

Ss ~ I dT I * X ~ fa"? + 2^ F ^ + tcF ^) ' 
where e 2 ~ m, a < 4, b and c are constants of order 0(1). F^ u = d^A s u — d u A s ^ and 9/2 is 
the "uniform" phase field for the Schwinger bosons. Notice that the gauge symmetry of the 
'uniform' gauge field A is broken by < Ajj >= ^ in SBMFT. The existence of stable 
vortex solution in SBMFT is tied with the A field as in usual superconductors. 

To derive vortex dynamics we concentrate at S u and introduce vortex 3-current j v in 
the phase field, 3^9 — > 3^9 + 3^9', where 9' is multivalued and j v = e^xdydxO' ^ 0. In the 
London limit where we treat <fi as constant a duality transformation can be performed where 
we can integrate out the 9 field to obtain an effective action S v for vortices fl4fl , 



S v = J dr J d 2 x x a) 2 + ia.(f - V x A)j + 



(4a) 



where V x a ~ (V6' — A) and corresponds to the transverse part of supercurrent field in 
superconductors. Notice that the amplitude field goes to zero at the center of vortex. As 
a result there exists an additional contribution to the vortex action, ~ {(energy needed 
to create vortex core, e v )x (length of vortex trajectory in space-time)}. For N vortices, 

~e v j:fdk = e v j:j d Jl + -L&*, (4b) 

i=l J i=l J V c v ar 

where Xj(r) represents the trajectory of the ith vortex in Euclidean space-time, c v is the 
limiting velocity. Notice that the suppression of <fi field at vortex core also couples vortices 
to the staggered amplitude and phase fluctuations through the term (— ^)q 2 in S s f9j. We 
shall first consider S v in the following. 

The dynamics of a single vortex can be obtained by minimizing at real time. We 



find that describes relativistic particles with energy E = ^ye v , where 7 = lA/1 — % and 

V c v 

v is the vortex velocity. In the absence of V x A term the particles have 'charge' (vorticity) 
and interact with each other through an effective U(l) gauge field a. For boson vortices 
the corresponding quantum field theory is a relativistic theory of scalar electrodynamics 
with charged bosons (vortices) ||14|| . In the presence of trapped magnetic flux inside vortex 
core (V x A — j v ~ 0), the electric 3-current j v is screened and the bosons decoupled from 
the gauge field a. The resulting theory is a relativistic theory of bosons with short-ranged 
interactions, as is in the case of vortices in usual superconductors |14] . 



To derive dynamics for the topological spin excitations we assume that once the boson 
spins are bound to the vortex core, their spatial degree of freedom is quenched and the 
only modifications to the pure vortex action (|J) are: (l)the vortices now carry spin indices 
m = —S, —S + 1, S, and there are 2S + 1 spin component of vortices, and (2)vortices 
become fermions when 2S is odd. In particular, since the vortex action (|j) is Lorentz 
invariant, the dynamics of the topological spin excitations must be described by (2 + l)d- 
relativistic field theories of bosons when S is integer and (2 + l)d- relativistic field theories of 



fermions when S is half-integer. To proceed further we examine more carefully the symmetry 
constraints in SBMFT. The existence of two lattice sites per unit cell in SBMFT implies that 
there are also two quantum fields ip^ and ip B in the continuum theory, representing spin-cr 
vortices centered at A- and B- sublattice sites, respectively Q. Under reflection or rotation 
by 7r/2 around center of a square plaquette, the A- and B- sublattices are interchanged and 
correspondingly also ip^ and ip B . Notice that we have considered finite temperature where 
< Z >= in our discussion and correspondingly parity (space-time reflection) symmetry is 
unbroken. In order to describe coupling of the bound boson spins at the vortex center to 
the fictitious gauge fields A and A s , the quantum fields ip£ and ip B must also carry 'charge' 
and are complex. 

With these kinematics constraints, the quantum field theories for the topological spin 



excitations can be determined |15[]. In the case of integer spins where the topological spin 
excitations are bosons, we construct quantum fields = ipa =t ''Pa which are eigenstates 

of Ti / 2 rotation and reflection with eigenvalues ±1. The dynamics of the tp^~^ fields are 
described separately by relativistic theories of complex scalar-fields. In the case of half- 
integer spins where the spin excitations are fermions the corresponding theory which respects 
parity is a theory of Dirac fermions in (2 + l)d ||15| . In the following we shall concentrate at 
the case of half-integer spin systems. 

In this case, the topological spin excitations are Dirac spinors described by four- 
component spinor field with 



MB) i 



\ 



(5) 



where ip^ B ^s are two component fermion fields needed to describe positive and negative 
energy solutions of the Dirac equation. In terms of ip a the effective Lagrangian which 
transforms correctly under parity is 



(6) 
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where 7^'s are usual 4x4 Dirac matrices in (2 + l)d with fj, = 0, 1,2. A is the 'uniform' 
gauge field. Notice that the staggered gauge field A s decouples from the ip a field in this level, 
as can be checked easily from transformation of L e ff under staggered gauge transformation 
^e id , and ip B ~^ ip„e~ l6 . Notice also that the spin degrees of freedom a appear as 
internal degrees of freedom for the Dirac fermions. 

It is interesting to compare the effective action L e // with the flux-phase |T]J of the 
fermionic mean-field theory of Heisenberg model in the case S = 1/2. In both cases pockets 
of Dirac fermions describe the fermionic spin excitations at low energy. The number of 
species of Dirac fermions are the same in both cases - there are four 'half-pockets' of Dirac- 
fermions around (k x ,k y ) = (±^,±^) in the flux-phase, and there are two 'full-pockets' 
of Dirac fermions coming from two independent combinations of ij)^ and ip^ in the present 
theory. The position of the 'Dirac-pockets' in /c-space cannot be determined with certainty in 
our effective theory. Nevertheless, to make comparison with the flux-phase we shall assume 
that the fermion pockets in our theory are centered around (k x , k y ) = (±| , ±| ). With this 
assumption the only difference between L e ff and the flux-phase is that there is a gap e v in the 
Dirac-fermion spectrum in our theory. In fact, in the limit e v = it can be checked directly 
that L e ff describes the continuum theory for the flux phase. The strong similarity between 
our effective theory and flux-phase leads us to speculate that the flux-phase in fact describes 
a new spin-disordered phase of 2D Heisenberg antiferromagnets where antiferromagnetism 
is destroyed by driving the fermion mass-gap a v to zero. We shall discuss this scenario in a 
separate paper. 

Next we consider the experimental consequences of the topological spin excitations. We 
shall first consider S = 1/2 as example. In terms of ip, The spins carry by the Dirac fermions 
in the S = 1/2 case are described by the operators 

S A{B \x) = \ £ £ (: ^ B) {x)B mm 4iS\x) :) (7) 

m,m' k=l,2 

where S A (x) and S B (x) are spin operators in the A- and B- sublattices, respectively and 
a are Pauli matrices, m,m' = ±|. : O := O— < O >g are normal-ordered operators. 
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The 'uniform' and 'staggered' spin operators can be constructed from S A ( B \ where m{n) = 
S A + (—)S B describes spin fluctuations around momenta q = (0,0) and q = (vr, 7r). Notice 
that because of the constraint that there is one spin per site in the original Heisenberg 
model, the Dirac fermions together with the original boson spins must satisfy the constraint 

e((E : cwcw o + z:(x)z:(x)) = 1, (8) 

tr \ fc=l,2 / 

which introduces additional coupling between boson spins and topological spin excitations 
not included in L e ff. 

Experimentally the Dirac fermion spectrum can be observed directly in neutron scatter- 
ing experiment at energy uj > 2e v . Their contribution to dynamic structure factor S(q,ui) 
at q ~ (0,0) and q ~ (n, n) can be calculated directly using Eqs. @ and (0). We find 
that their contributions are similar to S(q, uS) calculated from flux phase, except that a gap 
~ 2e v is found in the spectral function. The existence of Dirac fermion spectrum also af- 
fects spin correlation at lower temperatures. To see that we consider the antiferromagnetic 
spin-correlation length £(T). In pure SBMFT, £(T) is given at low temperature by |I(J 



aT) = V2^exp(^-n B ), (9) 

where n B is the density of bose-condensed spins in SBMFT at T = 0. In the presence 
of Dirac fermion spectrum, the system gets more disordered at T ^ because of thermal 
effect associated with the extra degrees of freedom and £(T) decreases. The effect can be 
estimated in a mean-field approximation using the constraint (§) where the average number 
of Schwinger boson per site < n SB > is given by 

< n SB >= 1 - J2 <: ^ka : > • 

fc=l,2,<r 

It is straightforward to show that at low temperature, < usb >~ 1 — 2exp(— and 
the leading correction at low temperature to £(T) is obtained by replacing ris — > ns — 
2exp(— jr) m Eq. (^). The reduction in ^(T) continues until T ~ e v , when the number of 
thermally excited fermion spin excitations becomes large, and the spin correlation becomes 
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qualitatively different from the prediction of SBMFT or NLoM. For general spin value S, 
e v ~ JS 2 and the temperature at which the NLaM description becomes invalid occurs at T ~ 
JS 2 which is deep inside the 'quantum critical' regime for large value of S. For small value of 
S where there is no clear separation between energy scales SJ and S 2 J, the whole 'quantum 
critical' regime may be washed away by the existence of topological spin excitations. Such 
scenario was indeed seemed to be observed in Monte Carlo simulations of S = 1/2 2D 
Heisenberg antiferromagnet f§-g| where the 'quantum critical' regime predicted by NLcrM 
description seems to be missing. Notice that for integer spin systems the contribution to 
dynamics structure factor S(q,u) from topological spin excitations are quite different from 
half-integer spin systems because of different statistics. However their effect on £(T) should 
be qualitatively similar. 

Lastly we discuss the effects of gauge field A. In the T — > limit where bose-condensation 
of boson spin takes place (< Z >^ 0) the gauge field A becomes confining |L7j and the 



topological spin excitations are confined in pairs by linear confining potential. At finite 
temperatures the confining potential is effective up to length scale ~ correlation length £ (T) 
and the effect of confinement is expected to be strong at low temperature when the system 
is at the 'renormalized classical' regime. As a result the topological spin excitation spectrum 
at q ~ (0, 0) is expected to be strongly modified from the free Dirac fermion prediction (for 
half-integer spins) at low temperature. However the behaviour of topological spin excitation 
spectrum at q ~ (jr, it) is determined by short distance behaviours of spin pairs and should 
not be affected strongly by confinement. As a result, we expect that our predictions for 
S(q,u) at q ~ (tt, 7r) and £(T) remains qualitatively valid. 

Summarizing, using SBMFT, we show in this paper that stable topological spin excita- 
tions exist in Heisenberg antiferromagnets at two dimensions. The topological spin excita- 
tions are described by relativistic quantum field theories of complex scalars when spin value 
S is an integer, and are Dirac-fermions when S is half-integer. Theoretical and experimen- 
tal consequences of these excitations are discussed where we point out that the existence 
of these spin excitations may be the reason for the anomalous results observed in recent 
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experiments || and Monte-Carlo simulations In particular, these excitations can be 

observed directly in neutron scattering experiments, and would provide a direct test of our 
theory. 
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